Lower energy bounds for quantum lattice 
Hamiltonians 

Dean Lee^and Nathan Salwen^ 



1 2 
deanjee@ncsu.edu, nathan_salwen@ncsu.edu 



Department of Physics, North Carolina State University, Raleigh, NC 27695 
Lcsu.edu, ^nathan_salwe 

February 1, 2008 



Abstract 

We derive general lower energy bounds for the ground state energy of 
any translationally invariant quantum lattice Hamiltonian. The bounds 
are given by the ground state energy of renormalized Hamiltonians on 
finite clusters. 



1 Introduction 

When computing the low energy eigenstates for some quantum Hamiltonian, H , 
one usually extracts an energy estimate by calculating the expectation of H for 
some approximate wavefunction. In the case of the ground state this estimate 
also provides an upper bound on the energy. The calculation is usually repeated 
with systematically improved wavefunction approximations and the final result 
is determined by extrapolation. While in principle this procedure should work, 
in practice computational difficulties can block improvement beyond some point. 
One common culprit is the fcrmion sign problem. Although the sign problem 
can be delayed with some success, as shown in [Q, in the end it will eventually 
lead to exponential scaling of computational effort. In such cases it is quite 
useful to have some idea of how far a given estimate is from the exact answer. 

Since in most cases the result in hand provides an upper bound to the ground 
state energy, what is needed is a reasonably accurate lower bound. In this letter 
we state and derive a general lower bound for the ground state energy of any 
translationally invariant lattice Hamiltonian. Our result is very general and 
applies to any number of dimensions, size of lattice, type of particle, or form of 
local interactions. In particular this includes single-site terms, nearest and next- 
nearest neighbor hoppings, and gauge plaquettes. The only restriction is that 
the interactions are localized and do not extend throughout the entire lattice. 
The lower bound is stated in terms of the ground state energy for a renormalized 
Hamiltonian defined on a cluster subset. As we find in the example described 
below, the error appears to scale as the ratio of the boundary size to the bulk 
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size of the cluster. Since the lower bound is defined on clusters, it is most 
convenient for position-space computational schemes which already partition 
the lattice into clusters. This includes density matrix renormalization group 
methods and Monte Carlo simulations using Hamiltonian partitioning 

Our letter is organized as follows. We begin by proving a general inequality 
involving two Hamiltonians which are equivalent within a projection subspace. 
We then apply the general inequality to a simple but illustrative example in- 
volving a one-dimensional periodic chain Hamiltonian and a non-periodic chain 
Hamiltonian. We then state and prove our main result, the lower energy bound. 
To conclude we check the lower bound numerically for our one-dimensional ex- 
ample and note the scaling of the error with cluster size. 

After posting this preprint to the archive, it has come to our attention that 
the lower bound techniques presented here have already been discussed in the 
literature /@/-^. 

2 General inequality 

We prove a simple but useful result from which the lower energy bound follows 
as a special case. Let H he a self-adjoint Hamiltonian. Let P be a projection 
operator that projects onto a subspace Vp. We assume that P commutes with 
H. Let H' be any self-adjoint operator whose spectrum is bounded below and 
satisfies 

PHP = PH'P. (1) 
Let E'q be the ground state energy of H' . For any vector v we find 

(Fw, HPv) = {v., PHPv) = {v, PH'Pv) = {Pv, H'Pv) . (2) 

Since 

(Pw, H'Pv) > E'„ {Pv, Pv) , (3) 

we conclude that 

{Pv, HPv) > E'q {Pv, Pv) . (4) 
Therefore the spectrum of H restricted to Vp is bounded below by E'q. 



3 Application to lattice Hamiltonians 

We now apply the general result (|4|) to a translationally invariant lattice Hamil- 
tonian. For clarity we consider a simple but sufficiently general example. Our 
example is a length N one-dimensional periodic chain with one type of fermion. 
We define Oi as a hopping term to and from the origin 

Oi = ajao + ajai, (5) 
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and O2 as the number operator at the origin, 

O2 = ajao- (6) 
Our Hamiltonian describes a free one-dimensional fermion, 

H ^ ci ^ aj^+ifln + a,\a„+i +0-2 ^ al^Un- (7) 

Ti=0,...,A'-l n=Q,...,N-l 

Let us now take a four-site non-periodic chain and the renormahzed Hamil- 
tonian H' , 

H' = f TnOi+f ^«^2. (8) 

n=0,l,2 n=0,l,2,3 

The reason for the renormalization of coefficients will be apparent shortly. We 
define Eq as the ground state energy of H' and take P as the projection operator 

n=0,...,N-l 

We note that 

PH'P = f ■ 3 • POiP + f • 4 ■ PO2P (10) 
= —PHP 

From (U) we conclude that the spectrum of H in the translationally invariant 
subspace is bounded below hy N ■ E'q. 



4 Lower energy bound 



We now state the general result, a lower energy bound for finite periodic lattices 
of any size or number of dimensions. Let be the operator corresponding 
with translation by the vector T^. We observe that a local operator such as 
a-s* transforms as T^a-s*T^^ — a^+v- Let G be any subgroup of lattice 
translations and \G\ be the number of elements of G. Let H he a. lattice 
Hamiltonian defined by 



H = 



E 

T^eG 



E CjT^OjT. 



(11) 



where the Cj's are coefiicients and Oj's are general quantum operators. Let H' 
be a modified Hamiltonian defined on a subset of lattice points. 



^ E ^ 

k—l,...,nj 



(12) 



The rij's are arbitrary positive integers and each ~T^j^k is a lattice vector (not 
necessarily distinct). Let E'^ be the ground state energy of H' . Since H' and 
are equivalent within the G-invariant subspace, we conclude from ^ that 
the energy spectrum of H is bounded below by \G\ ■ Eq. 
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5 Numerical example 



We now numerically check the lower energy bound for our one-dimensional 
chains. We use the notation Eq {01,02) for the ground state energy of the 
periodic chain Hamiltonian in (Q). This is trivially calculated in our example 
by momentum decomposition, 

E^{cuC2)^ J2 min(2cicos(2^) + c2,0) . (13) 

k=Q,...,N-l 

It easy to see that the lowest point of the spectrum in the translationally 
invariant subspace is also E^ {ci, C2). We take E'^^^jj^ , |j) to be the ground 
state energy of the non-periodic chain Hamiltonian of length M. We have 
renormalized the coefficients as prescribed in (|2|). In Tables 1 and 2 we show 
Eq (01,02) and the lower energy bound N ■ E'q '^{ ^^_^ , H) for several values of 
N, M, ci, and 02- The values for E'q^'{j§^, ||) were calculated by Lanczos 
diagonalization. 



Table 1. Lower bounds and energies for ci = — 1, C2 = 



N 


N-E'o" 


N-E'^ 


N ■ E'q^'^ 


N ■ E'a^"^ 


^0 


16 


-11.93 


-10.88 


-10.61 


-10.49 


-10.05 


32 


"23.85 


-21.75 


-21.23 


-20.99 


-20.31 


64 


-47.70 


-43.51 


-42.45 


-41.98 


-40.71 


128 


-95.41 


-87.02 


-84.90 


-83.95 


-81.47 



Table 2. Lower bounds and energies for ci = — 1, C2 = 1 



TV 


N -E'^^ 


N-E[f 


N ■ E'^^'^ 


N ■ E'^*^ 


^0 


16 


-4.63 


-4.08 


-3.90 


-3.76 


-3.53 


32 


-9.26 


-8.17 


-7.79 


-7.52 


-7.00 


64 


-18.52 


-16.33 


-15.59 


-15.05 


-13.96 


128 


-37.04 


-32.67 


-31.18 


-30.10 


-27.91 



We note that the error between the lower bound estimate and the actual energy 
scales as M~^, the ratio of the boundary size to the bulk size of the cluster. 
This is also expected in the general case. 

6 Summary 

We have derived lower energy bounds for any translationally invariant lattice 
Hamiltonian. Our result is general and applies to virtually any local quantum 
lattice Hamiltonian relevant to particle, nuclear, or many-body physics. The 
lower bound is given in terms of the ground state energy for a renormalized 
Hamiltonian defined on a subset of lattice points. We find that the error scales 
as the ratio of the boundary size to the bulk size of the cluster. Given the 
simplicity, generality, and convenience of the result, we hope this lower bound 
estimate will be useful in many future computational lattice calculations. 
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